Quantum simulation of expanding space-time with tunnel-coupled condensates 
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We consider two weakly interacting quasi- ID condensates of cold bosonic atoms. It turns out that 
a time-dependent variation of the tunnel-coupling between those condensates is equivalent with the 
spatial expansion of a one-dimensional toy-Universe with regard to the dynamics of the relative 
phase field. The dynamics of this field is governed by the quantum sine-Gordon equation. Thus, 
this analogy could be used to 'quantum simulate' the dynamics of a scalar, interacting quantum 
field on an expanding background. We discuss, how to observe the freezing out of quantum fluctu- 
ations during an accelerating expansion in a possible experiment. We also discuss an experimental 
protocol to study the formation of sine-Gordon breathers in the relative phase field out of quantum 
fluctuations. 
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Figure 1. (Color online) Quantum simulator, (a) Sketch of 
the quantum simulator consisting of a pair of tunnel-coupled 
ID bosonic condensates. The tunnel amplitude t± is modu- 
lated time-dependently. It can be related to the scale param- 
eter a (b) describing the (exponential) expansion of a 1+1 
dimensional 'Universe'. The laboratory time (condensate co- 
ordinate x) is identified with the conformal time 77 (co-moving 
coordinate) (see main text). The relative phase field between 
the condensates (j) represents the scalar quantum field to be 
simulated, (c) Semiclassical evaluation of the physical en- 
ergy density of the simulated quantum field (here displayed 
in physical coordinates Xph = ax). One can observe the struc- 
ture formation out of zero-point quantum fluctuations during 
the expansion (see main text). 



The recent progress in coherently controlling systems 
of cold atoms (e.g., pW] , see also [5] and references 
therein) , stimulated a lot of research concerned with em- 
ploying these experimental systems to 'quantum simu- 
late' prototypical quantum many-body models (see e.g., 
[6l[7]). Particularly fascinating are ideas concerned with 
simulating quantum many-body physics on curved space- 
time ('analog gravity') (see, e.g., [8:-il4j and [15j, as well 
as references therein) connecting concepts and techniques 
from cosmology and condensed matter (see [16j). 

Here, it is argued that a pair of tunnel-coupled, quasi- 
ID, bosonic condensates (as realized in [17J and related 
experiments) can be employed for simulating an interact- 
ing, scalar quantum field on top of an expanding 1+1 di- 
mensional space-time (Fig.[l^,c). This scalar field is rep- 
resented by the relative phase field between the conden- 
sates. As argued in [T8"/19j, at low energies, its dynam- 



ics is described by the quantum sine-Gordon model. It 
turns out that in this setup, one can simulate the expan- 
sion of the ID toy-Universe simply by varying the tunnel- 
amplitude according to a suitable protocol (Fig.[T]3). In 
the experiments (e.g., [17J) the tunnel-amplitude itself 
can be largely tuned and the field dynamics can be di- 
rectly visualized by means of matter-wave interferome- 
try. The 'quantumness' of the relative phase-field dy- 
namics depends on the interaction strength within each 
condensate (which is easily tunable, e.g., by adapting the 
ID condensate density). As far as we are dealing with 
non-linear dynamics, we consider the weakly interacting, 
'semiclassical' limit (e.g., [T91j2l] ) and employ the trun- 
cated Wigner approximation (TWA) (e.g., L2QJ). In con- 
trast, the proposed quantum simulator would allow to 
explore also the deep quantum regime - in fact, this is its 
main purpose. 

In the following, we start with deriving an effective 
Hamiltonian description of an interacting quantum field 
on an expanding 1+1 dimensional space time. In a second 
step, we introduce the experimental 'quantum simulator' 
- the tunnel-coupled condensates - and eventually discuss 
two of the possible effects, which could be explored. 

The first effect is the well known 'freezing' of quantum 
fluctuations and the related 'cosmological particle pro- 
duction' during an accelerating expansion. This purely 
linear, though most fundamental effect is made responsi- 
ble for the structure formation in the very early universe 
[22J. In the present experiments, one can only perform 
single measurements of the relative phase field per run. 
Fortunately, this mode freezing also manifests itself in 
the spatial fluctuation spectrum of the field and thus, 
in principle, is detectable with single measurements per 
run. Although, strictly speaking, there is no need for 
a quantum simulation of the exactly solvable linear dy- 
namics, observing the freeze-out of quantum fluctuations 
in the experiment would nevertheless be exciting (see for 
instance ilOj) and constitute an important check on 
the setup. 

The second feature is the generation of localized, 
macroscopic structures during the expansion out of quan- 
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turn fluctuations. This pattern formation involves the full 
non-linearity of the underlying sine- Gordon fleld theory 
and was also observed in the static case [19j. In contrast 
to [19], for an exponentially fast expansion this happens 
only for small enough expansion rates. At large expan- 
sion times, these patterns seem to turn into standing 
sine- Gordon breathers simply drifting apart from each 
other. We argue how to detect signatures of this 'Hub- 
ble' drift experimentally. In cosmology, e.g., dealing with 
the preheating following inflation, excitations like this 
(e.g., in the scalar inflaton fleld) are sometimes denoted 
as 'oscillons' (e.g., [23H28] ). A full experimental quantum 
simulation would allow investigating the formation and 
persistence of these excitations on an expanding back- 
ground, even in regimes where quantum effects become 
very important for the dynamics. 

Quantum field on curved space-time. - The space-time 
action of a classical, scalar fleld theory in 1+1 dimensions 
is given by [22] (c = 1, = 1) 

S=ll dx^^\detg^,\[id^x)9'"'id,x) - ^V{x)] , (1) 

where g^jy denotes the metric (with g^vd^^ — ^tJ^X = 
dx/dx^ and V{x) is an arbitrary potential. We are in- 
terested in an homogeneous, spatially expanding space- 
time described by the Friedmann- Robertson- Walker met- 
ric (FRW) (see for instance [22j, ds'^ = g^ydx^dx^) 



ds' 



dr'^ — a^{r)dx'^ . 



(2) 



Here, x denotes co- moving coordinates which are related 
to physical coordinates Xph via the scale parameter a(r) 
as 2Jph ~ dx . The time r denotes the cosmological time, 
i.e., the proper time of a co- moving observer. In reality, 
the dynamics of a is determined by Einstein's equations. 
Here, we treat a(r) as a given function of time, which 
will be specifled below. 

In 1+1 dimensions, the Lagrangian corresponding to 
the action Eq. ([T]) takes an intriguingly simple form using 
the so-called conformal time r] with dr] = dr/a: 



dx 



(5,x)'-axf-2a2(r,)y(x) . (3) 



Note that light geodesies {ds^ = 0) are described by 
x{r]) = ±r] + const. The quantization of the fleld theory 
follows the standard prescription (e.g., [22j). First, from 
Eq. ([3|, we construct the canonically conjugated fleld 
H^(x, 1]) = = drjX' ^ second step, we promote the 

flelds to operators demanding that [x(^, ^x(^^ ^)] = 
i6{x — x^). Eventually, we can switch to the Hamiltonian 
formulation introducing the time-dependent Hamiltonian 



H{v) 



dx 



tll + {d,x) +2a\v)V{x) ■ (4) 



Note that all effects of the expanding space-time are now 
encoded in the time-dependence of (77) (cf. [T6|f29[ l3Q ] ). 



Tunnel coupled condensates as quantum simulator. - 
Here, we propose a quantum simulation of the fleld 
x{x^r]) for the special case of a sine-Gordon potential 
V = — ttiq/^"^ cos(/3x). This potential has several inter- 
esting properties: The corresponding fleld theory is inter- 
acting and integrable. Second, the sine-Gordon potential 
appears in the so-called 'natural inflation' scenario [31j. 
Third, the sine-Gordon potential supports the formation 
of 'quasibreathers' [19j (in the cosmology literature de- 
noted as 'oscillons', e.g., [25] |26l [28l |32j and references 
therein). 

The quantum simulator (Fig. [l^) consists of two 
tunnel-coupled quasi- ID condensates of cold, bosonic 
atoms (e.g., [17J) with a time-dependent tunnel ampli- 
tude t±. The laboratory time is identifled with the con- 
formal time 77. At low energies, the dynamics of the rel- 
ative phase fleld l3(j){x^t)/^/2 = (^1 — ^2)/\/2 can be 
described by the quantum sine-Gordon model [18j (the 
sound velocity Vs = I) 

rL 



dx 



2rrr{ri) 



cos 



.(5) 



The relative phase fleld and the relative density varia- 
tions a/2/3~^II = ^^^^ form a canonical pair. The tun- 
nel amplitude enters the mass term m?{r]) = 2^'^ potj_{r]), 
where po is the mean density per condensate. The Lut- 
tinger liquid description should be reliable as long as the 
typical lengthscale of Eq. ([5|, set by \^hvs/rri, is much 
larger than the healing length of the condensates 
[m IE]- This can always achieved by choosing a suffi- 
ciently small tunnel amplitude t±{r]). The parameter /3 
is related to the Luttinger parameter K as P = ^J2i\ jK. 
For weak interactions /3 <C 1. This is the limit we are 
considering here. It can be shown that /3 plays the role 
of Plank's constant ^20j and /3 <C 1 corresponds to the 
semiclassical limit of the quantum sine-Gordon model 
(see also, e.g., [33j). The analysis here (as far as we are 
dealing with non-linear dynamics) is based on the semi- 
classical TWA. However, in the experiment one can go 
deep into the quantum regime corresponding to larger /3 
(a rather broad range of values up to K ~ 50 is real- 
izable, e.g. [34J). The following identiflcations connect 
the quantum simulator and the quantum fleld theory on 
an expanding background. Identifying the flelds x ^ (\) 
and rv?{r\) = mQa'^{r]) [and thus a{r]) ^ \/tjJj]Y/t±]fi)], 
the dynamics of the relative phase fleld simulates x 
conformal time and co-moving coordinates. In the re- 
mainder, we will always argue in terms of the fleld ^, 
i.e., we analyze Eq. ([5|. 

Scale parameter and initial state. - We consider an 
exponential expansion a(r) = e^^ with the 'Hubble 
constant' H = d{r)/a{r). Choosing a(0) = 1, the 
conformal time is given by 7^ = [l — e~^^] and 

correspondingly a{r]) = [1 — Hr]]~^ (see Fig.[l]3) with 
T] G [0,7^/] and rjf < 1/H. At 7^/, the expansion ends and 
a{r] > rjf) = a/. The dimensionless parameter H/mo 
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Figure 2. (Color online) Mode 'freezing', (a) Plot of the 
spectrum at a = 19.6, 166.8, 1000 (blue lines, top to bottom). 
Here, H/mo = 2.6 and mo = 0.05. One observes the spec- 
tral tilt (dashed, blue line) resulting from the freezing out of 
modes. Dashed lines show the corresponding ground state 
spectra for a = 1 and a — af. (b) 'Freezing out' of modes 
(here k/afmo = 0.1). While the suppression of the ID field 
fluctuations at later times is non- universal (e.g., it depends 
on the spatial dimension), the mechanism of the expansion- 
induced mode freezing is the same as in higher dimensions, 
(c) Sketch of the horizon crossing 



corapares the expansion time-scale and the typical inter- 
nal time scale (at short times) of the system and plays a 
crucial role throughout the following. 

At 7^ = 0, we start in the ground state of massive 
phonons with a small mass mo- In particular, mo is 
chosen much smaller than the UV cutoff 1 ~ 1/^/^ (all 
TWA simulations are performed on a lattice with lat- 
tice constant set to one, while keeping mQa{r]) <C 1 
throughout the whole simulation). For finite system 
size and /3 <C 1, the center-of-mass mode (COM mode) 
<l>(7^) = dx ^(x, rf) needs some special attention. For 
L/ih ^ we can treat the COM mode classically such 
that (p^x^rf) ^{v) -\- S(j){x^r]). Before, the expansion 
starts, the COM mode is tuned to some value ^(0) ~ ^0 
with e [0,7r] {^'{0)/mo c^O). As argued in [19j, such 
an initial state can be achieved by slowly splitting a sin- 
gle condensate followed by applying a potential gradient 
between the condensates to tune <l>o. 

Freezing out of quantum fluctuations. - One of the fas- 
cinating results of modern cosmology is that the struc- 
ture formation in the very early Universe seems to have 
been seeded by quantum fluctuations [22j. This result is 
truly amazing, as on cosmological scales zero-point fluc- 
tuations are tiny. However, it seems that an exponential 
expansion of the very early Universe (inflationary stage) 
led to a 'freezing' of quantum fluctuations and stretched 
them to cosmological scales. One can reformulate this 
basic mechanism in a condensed matter language (e.g., 
|16|). In this terminology, the inflationary expansion cor- 
responds to a rapid, non-adiabatic 'quench' (see, e.g., 
[35J) producing a large number of excitations ('cosmo- 
logical particle creation', cf. |10|). 

A very similar effect should be observable in the con- 
sidered 1+1 dimensional toy-Universe. For this pur- 



pose, we consider the case ^0 — 0- For small enough /3 
and finite system size, one can safely expand the cosine- 
potential to lowest order m^(77)/3~^(l— cos /3(^) ^ ^ 
yielding a theory of massive phonons (cf. [35H38] ) 
on an expanding 1+1 FRW space- time. According to 
Eq. (|4|, the dynamics of the modes 4>k [4>{^^v) = 4* + 
Efe/o e'^'^Mv)] follows as 

^'ar])^[k'^my{ri)]Mr])^0. (6) 

The general solution of Eq. (|6| is given in the Sup- 
plement. For \k\/aH — \\ (where u = 
^^/Y^^^moJlT^)^ the modes evolve freely. However, for 
\k\/a{r])H ^ - \ \ (cf. Fig.|2]D,c), the modes 'freeze'. 
Most importantly, this 'freezing' also manifests in the 
spectrum {(j)k{'n)^-k{v)) (cf- [9. 16,^). This, in prin- 
ciple allows to observe this effect with a single measure- 
ment of the relative phase field ^(x, rf) per run. For the 
considered protocol (restricting to 77 < r^j), we obtain for 
all modes with \k\ ^ H\v^ — ||,mo (starting inside the 
"horizon") and \k\/a{r]) <C H (see Fig.|2}i) 

1 H^^~^ 

{Mv)4>-k{v)) oc ^13^^^. (7) 

Modes with \k\/a{r]) ^ Hlv^ — || remain close to the 
initial ground state of massive phonons with mo yield- 
ing {(l)k{r])(j)-k{'r])) — 1/2|A:|. Here, we consider a fast 
expansion with H/mo > 2. In the limit H/mo ^ co, 
u approaches 1/2 corresponding to a 'sudden quench', 
which leaves the spectrum unchanged. In a possible ex- 
periment, one could detect the different power-laws by 
probing the longitudinal (relative) phase coherence ^39j 
of the condensates on different length-scales (very simi- 
lar to [40j). It can be easily checked that the finite a'{0) 
does not influence the 'particle production' (indicated by 
the deviations from the instantaneous ground state spec- 
trum [see Eq. ([t])]), which happens happens during the 
evolution. 

Formation of breathers out of quantum fluctuations. 
- While the freezing out of quantum fluctuations is a 
purely linear effect, we now discuss a feature, which heav- 
ily relies on interactions. From now on, we consider a 
slow expansion H/mo 1 and finite < <l>o 7r//3 
(see Fig. [3^). At short times 77 > 0, the global rel- 
ative phase ^{t]) performs Josephson oscillations (see, 
e.g., [19j). However, it is well known that these are 
parametrically unstable against small, spatial fluctua- 
tions (e.g. (El im |42j ) . Linearizing around ^, one finds 
^k + [^^ + m^a^ {t]) cos /3^]4)k — 0. In the static case 
[19j, this parametric drive leads to ^^kiv) ^ e^ik,^o)r] 

with r J^^sin^(/3<l>o/2) - P/m^ f42] (displayed for 

/3^o ^1)- In turn the fluctuations (Scj)'^) grow and at 
some point the linearization breaks down. In the semi- 
classical limit /3 <C 1, it was demonstrated [19j that the 
non-linearity of the sine- Gordon equations leads to the 
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Figure 3. (Color online) Emergence of breathers, (a) Plot 
of a single TWA run in co-moving coordinates and conformal 
time {/3 = 0.1, /S^o/tv = 0.65, H/mo = 0.02). One observes 
the damped Josephson oscillations of the COM mode and 
the emergence of breather-like excitations out of parametri- 
cally amplified phonons. As the COM mode is damped away 
due to the expansion, standing sine-Gordon breathers are re- 
leased. These can be tracked numerically (green dots), (b) 
Time-evolution of these breathers tracked at Sr] = denot- 
ing the time of tracking (77 > 77/) as a function of QbSt] (Qb 
is the numerically obtained breather frequency). The spatial 
center of the breathers is shown. Color code: Distribution 
function; Blue, solid line shows the temporal evolution of a 
standing breather (corresponding to the mean breather fre- 
quency) . (c) Density of created breathers {J\f is the number of 
tracked breathers) obtained by tracking breathers at a/ = 10. 
Here, L = 1000, mo = 0.075 and /3 = 0.1. The white, dashed 

line indicates where /3max^ (^0^) = 0-2 (see main text). 



formation of localized patterns in the field ^(x, 77). These 
patterns were identified with 'quasibreather'-solutions of 
the classical sine- Gordon equation [19], which in con- 
trast to usual breathers have a finite lifetime. The 
'quasibreather'-solutions could also explain the formation 
of these excitations out of phononic quantum fluctua- 
tions. 

Here, for a slow, exponential expansion (in proper time 
r), one observes the formation of similar excitations only 
for mo/H exceeding a certain, almost sharp threshold 
(Fig.jsj^, cf. also [28j), which depends on /3 and ^q. To see 
this, first note that as long as the linearization of the sine- 
Gordon equation applies, all field fluctuations are sup- 
pressed as ^ (including the COM mode). Fur- 
thermore, in the course of time, resonant modes can get 
shifted out of resonance as a consequence of the expan- 
sion. For a slow expansion, replacing k ^ k/a and ^0 
^o/Va (cf. [Hgg), one finds that {^k{r])^-k{r])) 



, which is valid, as long as 



a exp 

ImT = 0. The non-linearity and thus the formation 
of localized patterns kicks in only if the fluctuations 




Figure 4. (Color online) Plot of the correlation function 
^00(0,77). Clearly, one observes the exponentially growing 
fluctuations (parametric resonance). At some point the non- 
linearity of the sine-Gordon model leads to the formation of 
breather- like structures (around the main peak, cf. Fig.[3^). 
While the COM mode is damped by the expansion, these 
structures turn into standing sine-Gordon breathers with a 
constant amplitude (cf. Fig.jsjD), which, in physical coordi- 
nates simple drift away from each other. This 'Hubble expan- 
sion' is mirrored by the linear decay of the correlation func- 
tion as oc (1 —77 if) at late times. Once the expansion stops at 
77/, the correlation function becomes almost constant. Here: 
13 0.1, H/rriQ 0.02, tuq 0.05, L 800 and /3$o O.Ttt. 
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exceed a certain value for some r] (see Fig. [sj: and 
[43J ) . Numerically, we find that this value is of the order 
(9(10-1). 

Close to the creation threshold, once created, these 
breather-like excitations persist at the position, where 
they were 'born' out of quantum fluctuations. This is 
in contrast to the 'quasibreathers' observed in the static 
case [19j (cf. also [32j). In the long-time limit, the ho- 
mogenous part of the field $ is damped away cx \l\fa. 
We find good numerical evidence that the localized ex- 
citations, however, are robust against the expansion and 
can be well described as standing (classical) sine-Gordon 
breathers 0b (see Fig.[3]3). Their typical distance is set 
by the maximally amplified wavelength before the non- 
linearity sets in, ending the parametric amplification. 
It seems that at late times (a ^ 1), the only effect 
of the 'adiabatic' expansion {H/ttiq <C 1) on breathers 
is a trivial shrinking of the breather period and width 
(both (X l/mQa{r])) in co-moving coordinates, while their 
amplitude stays approximately constant. This can be 
understood realizing that the amplitude of a classical 
sine- Gordon breather /30b is solely determined by the 
breather parameter G [0,7r/2] (max/30B = (27r — A^p)^ 
see e.g., [44J). In the quantum sine-Gordon model, this 
parameter gets quantized [45j, i.e., it is promoted to a 
quantum number. However, it is well known that for 
a slow change of system parameters (by slow, here, we 
understand 0.^^80.^/97] <C 1, where I^b = rnQa(r]) s\n Lp 
is the instantaneous breather frequency), quantum num- 
bers (and thus the breather amplitude) are approxima- 
tively preserved (cf. [46j). While for large a ^ 1 the 
number and amplitude of breathers remain constant (per 
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experimental run), they simply move apart from each 
other in physical coordinates. This 'Hubble expansion', 
e.g., can be observed in the (experimentally accessible) 
equal time correlation C00(x,7^) = (^(x, 7^)^(0, r^)) — 
(^(0,7^))^ (Fig.|4|. Under the assumption that at late 
times, the field (j) can be described as a set of independent 
(standing) sine-Gordon breathers with fixed amplitudes, 
one obtains that 



^0^(0,7?) (X 1 - r]H. 



(8) 



The linear suppression of is a direct consequence of 
the decreasing breather density in physical coordinates. 
From a condensed matter point of view the observation 
that a suitable protocol for the tunnel-amplitude pre- 
pares a state consisting of independent, standing sine- 
Gordon breathers is interesting by itself. While the anal- 
ysis here is based on semiclassical considerations (numeri- 
cally on the TWA) reliable for /3 <C 1, the proposed quan- 
tum simulator could for instance test the stability of clas- 
sical sine- Gordon breathers against quantum fiuct nations 
for larger /3 ~ 1 (cf. [47j). Furthermore, a quantum sim- 
ulation could give insight in the excitation of 'oscillonic' 
patterns (as discussed in the cosmology community, e.g., 
[26H281 I32]) in a scalar quantum field (such as the infiaton 
field or even the Higgs field) during a spatial expansion. 

Conclusions. - We demonstrate that tuning the 
tunnel-amplitude between a pair of tunnel-coupled ID 
condensates, the relative phase field can simulate an in- 
teracting quantum field on an expanding 1+1 space time. 
The proposed 'quantum simulator' should be realizable 
with present cold atom setups. As examples of the quan- 
tum many-body dynamics, which could be investigated, 
we discussed the freezing out of phonon modes and the 
creation of sine-Gordon breathers out of quantum fiuctu- 
ations during an exponential FRW-expansion. While the 
discussion here is restricted to the semiclassical limit of 
the underlying quantum sine-Gordon model, the 'quan- 
tum simulator' is meant to explore the deep quantum 
regime. 
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Supplemental material. - Here, we discuss the dynam- 
ics of massive phonons on the expanding 1+1 dimensional 
space-time in conformal time following Eq. (|6|. The clas- 
sical solution is given by 



(9) 



where u = ^ y^l — AttIq/H'^ and Jj^ denotes the Bessel 
function of the first kind. In the following, we restrict 
to H/mo > 2. The quantum mechanical solution fol- 
lows from Eq. ^ promoting the coefficients to oper- 
ators and by matching the initial conditions at = 
(a = 1), i.e., MO) = \ Ah + bl^) and 4>m = 



^^^^(6fe - bl^). Here, uJk{mo) = ^k^ + is the 

dispersion of massive phonons and 6/c|0) = (|0) denotes 
the ground state of massive phonon modes with /c 7^ 
and mass mo, in which the system is initialized). As 
Jiyix > {)) G M for liav = 0, we obtain 



Cl 



3—v /yj2uJk + iJ-j,^uJk/2 

j—yJv jvJ—y 



bk^albl,, (10) 



1 



/2ujkJ-v 



h^pfbl. 



(11) 



where we introduced the abbreviations j±jy 



L J 77=0 

tuation spectrum yields 



and J±j, = J±u 



. The fluc- 



(12) 



We are interested in modes starting within the horizon, 
i.e., with \k\ ^ H\u^ — l|,mo. Making use of the limits 



lim 



J.(|x|) = 
lim Ju{\x\) = 



2 ^TTiy 
sin(-— 

TTX ^ 2 

\x\' 



(13) 



one obtains for large times 77, such that \k\/a{r])H <C 1 
(and finite u): 



(14) 



